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A model of two-leg spin-S ladder with two additional frustrating diagonal exchange couplings Jd , 
J'o is studied within the framework of the nonlinear sigma model approach. The phase diagram has 
a rich structure and contains 2S gapless phase boundaries which split off the boundary to the fully 
saturated ferromagnetic phase when Jd and J'^ become different. For the S = \ case, the phase 
boundary is identified as separating two topologically distinct Haldane-type phases as discussed 
recently by Kim et al. ( cond-mat/9910023|) . 
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I. INTRODUCTION 

Low-dimensional spin models still continue to attract a 
considerable attention of researchers, both in theoretical 
and experimental aspects. Since the famous prediction 
of Haldane in 1983 [l] of different behaviour of Hciscn- 
berg spin chains with integer and half-integer value of 
spin S, which was based on a mapping to the nonlinear 
sigma model (NLSM), the NLSM approach was recog- 
nized as an important tool in studying spin systems and 
have found numerous applications (see, e.g., [g| for a re- 
view). Normally, the NLSM approach does not give good 
numerical results, however it is usually able to capture 
the topology of the phase diagram. 

During the recent upsurge of interest to spin lad- 
der models, several researchers have successfully applied 
NLSM to describe a A^-leg spm-S ladder |g-|]. In an 
essential similarity to the case of a single chain, it was 
found that for half-integer S, the ladders with even or odd 
number of legs N are respectively gapped and gapless. A 
natural question arose, namely, whether the properties 
of a gapped phase of, say, a two- leg spin-i ladder are in 
some sense equivalent to those of the Haldane phase of a 
spin-1 chain. 

Several arguments were given in favour of the positive 
answer to the above question. ^-g| Particularly, it was 
shown that by adding extra interactions one can intro- 
duce a suitable generalization of the pure ladder model, 
increasing the number of parameters in the phase space, 
and then one can find a path in this generalized phase 
space which smoothly (i.e., without crossing any phase 
boundaries) leads from the ladder model to a certain 
composite representation of a spin-1 chain [|6|,y|; more- 
over, it was demonstrated Mm that a two-leg S = ^ 
ladder has nonzero string order which is believed [^ to 
be a characteristic feature of the Haldane phase. 

On the other hand, it turned out that other general- 
izations may have very different properties; for instance, 
for the model of a "diagonal ladder" with additional 
equal-strength diagonal interactions, which also yields a 
composite-spin representation of a spin-1 chain [nOl one 



finds numerically that the "usual" ladder is separated 
from the composite-spin- 1 Haldane phase by a transition 
line |11^,&31 . An exactly solvable model exhibiting similar 
features was also constructed |l3|| . 

Recently, Kim et al. ]1J] have made an interesting ob- 
servation, noticing that there are actually at least two 
different definitions of the string order for a two-leg spin- 
i ladder (depending on whether one combines the S = ^ 
spins on the rungs or on the diagonals). Exploiting the 
analogy with the topological quantum numbers which 
can be introduced for short-range valence bond states 
on a square lattice |13|, they have conjectured that those 
two definitions of the string order distinguish between 
two different Haldane-type phases. This assumption was 
supported by the results of the bosonization study of two 
generalizations of the ladder model. 

In hope to get a better understanding of the physics 
of the spin ladder, and to search for possible new phase 
transitions, we find interesting to study the phase dia- 
gram of the generalized ladder model with unequal di- 
agonal couplings. We consider the model determined by 
the Hamiltonian 



H — Jl 2_^ /_^ Sa,iSa,i+l + Jr. 2_^ Si^iS2,i 

a— 1,2 i i 

+ 2_^^'^DSl^iS2,i+l + Jo^2,iSl^i+l) , (1) 

i 

where Sa,i are spin-S* operators at the i-th rung, a = 1, 2 
distinguishes the ladder legs. The model is schematically 
shown in Fig. nl At J^i = J^ = one recovers a reg- 
ular ladder, while at J^ = the model is equivalent to 
a zigzag chain with alternation of the nearest-neighbour 
interaction. Interchanging Ju and J^ is obviously equiv- 
alent to interchanging the legs of the ladder, so that it 
is sufficient to restrict ourselves to the Jd > J'jj case. 
The point Jd — J'd is in a certain sense special, since it 
allows an additional symmetry operation, namely, inter- 
changing the spins on every other rung is then equivalent 
to interchanging Jd and J^. 

The phase space of the model is three-dimensional 



and is determined by the three ratios of exchange con- 
stants, e.g., Jd/Jr, J'dI Jri Jl/ Jr- We will show that 
for Jd 7^ J'jj the phase diagram of the above model al- 
ways possesses 2S gapless phase planes, which split off 
the boundary to the fully saturated ferromagnetic phase 
at Jd = J'd- We consider in detail the most interesting 
case S — \ and show that the gapless plane is an exten- 
sion of one of the transition lines discussed in Ref . [Q . 
In the next section we briefly describe the mapping to 
NLSM, Sect. Ill contains the discussion of the results, 



and, finally. Sect. [V gives a brief summary. 



II. 



RESULTS OF THE MAPPING TO THE 
NONLINEAR SIGMA MODEL 



To map the model (nl) to a NLSM, we use the well- 
known technique of spin coherent states path integral; 
this technique is well described in reviews and textbooks 
0, and here we will not give a complete derivation but 
rather indicate only the main steps. We choose a four- 
spin plaquette as an elementary magnetic cell; then there 
are four classical ground states commensurate with this 
choice of cell, namely a ferromagnetic state (F) and three 
modulated states shown in Fig. Qand denoted as (A), (B) 
and (C). At n-th plaquette we introduce four variables 
in„, 1„, u„, v„, defined as the following linear combina- 
tions of the "classical" spin vectors (parameters of the 
coherent states): 



In — T7;(Sl,2ri-l + S2,2n-1 ~ Si_2n " S2,2n) : 
nin — T7;(Sl,2ri-l + S2,2n-1 + Si^2n + S2,2n) : 



Un == T7;(Sl,2n + S2,2n-1 " Si,2n-1 — S2,2n) • 
Vn = T7;(S2,2ri-l + S2,2n " Si^2n-1 ~ Si^2n) : 

which satisfy the following four constraints: 

m^ + 1^ + u^ + v^ = 1, 
(m + 1) • (u + v) = 0, 
(m - 1) • (u - v) = 0, 
(m + v) • (u -I- 1) = 0, 



(2) 



(3) 



Those variables we consider as smoothly varying func- 
tions of the space coordinate Xn = na when passing to 
the continuum limit; one should mention that the above 
ansatz is essentially similar to that used by Senechal [Q . 
The advantage of the ansatz (H) is that it conserves the 
total number of the degrees of freedom, which is impor- 
tant to avoid ambiguities in the mapping, as was recently 
realized on the example of inhomogeneous spin chains 

The order parameter for the four commensurate classi- 
cal ground state configurations F, A, B, C is respectively 



m, u, V, 1. Comparing the energies of those configu- 
rations, one may obtain a "draft" of the classical phase 
diagram which neglects presence of any incommensurate 
ground states; for the moment we are mainly interested in 
the commensurate antiferromagnetic part, and the con- 
ditions for the existence of spiral phases will be obtained 
later. One thus may treat m as a small fluctuation, and 
obtain different field descriptions starting from one of the 
configurations A, B, C. Massive degrees of freedom can 
be integrated out in a usual way, and in each case one 
obtains the final effective action in the form of a NLSM, 

A,s/h^j-J Jd^drlidtnf - {d^nf} 

+ ^JJd^dTn-{d.,nxdtn), (4) 

where n is the corresponding order parameter, and ^ = 
x/a, T = ct/a are dimensionless space-time variables, a 
being the lattice constant along the legs direction. For 
each of the classical "phases" A, B, C the coupling con- 
stant g and the topological angle 9 are given by the fol- 
lowing expressions: 
Phase A: Jr -h 2Jl > 0, J^ < Jr, J^ < 2Jl. 



9A 



Jr + 2 Jl 



6*^ = mod 27r, 



2S'VWa ' 

WA = \{JR + 2JL)[2JL-Ji 



(5) 



{Jd) 



k) 



4-'^ -^ - - Jr-.^ 

Phase B: Jr + J+ > 0, 2Jl < J^, '^Jl < Jr- 

JR + Jt „ 4^5 Jr 



>0. 



9B = 



D 



ys = 






2SVW^' " Jr 
Wb = 1{{Jr + JWd - ^Jl) - {Jof} 



(6) 



>0. 



Phase C: J? 



'D 



9c = 



Wc = J (4 



-2Jl>0, Jr<J^, Jr<2Jl. 
2Jl 



6c ^0 mod 27r, 



(7) 



2Jl){2Jl 



J 



+ {Jp? 1 

Jn ^ Jr 



>0. 



Here for the sake of convenience we have introduced the 
notations 

Jd — (Jd ± J'd) ■ 

The spin wave velocity for each case is given by c = 
2\/WSa/h. The inequalities define the boundaries of the 
domains of validity of the correspondent mapping (not 
all of them represent real phase boundaries, as will be 
discussed later). The boundaries defined by Wa,b,c = 
represent just the classical conditions for the transition 
into a spiral phase. One may observe that there is no 
spiral phase at J^ = 0. 



Phase F has to be considered separately, and it is easy 
to obtain its boundaries using the linear spin wave theory. 
There are two magnon branches with the energies 



£±(g) = -S{Jb. + Jd+ 2Jl) + 2SJLCosq 



1/2 



(8) 



and from the condition of positiveness of £± it is easy to 
obtain the boundaries of the F phase. They are deter- 
mined by the inequalities 



Jr + Ji > 0, Jr + 2Jl > 0, 



Wf 



2JL-Ji + {JD?/{JR + Jt>)>^- 



(9) 



At Wp = 0, £-{q) changes sign at once in a finite inter- 
val of wave vectors near g = 0, signaling the first-order 
transition, £-{q = it) vanishes at the line Jr + 2Jl — 0, 
and e+(g = 0) becomes zero at the line Jr + J^ = 0. 

One can see that only in the (B) case there is a non- 
trivial topological term, and the condition of gaplessness 
6 = {2n + l)7r yields 



Jd = 



2n + l 
45 



(Jfl + J+), 71 = 0,1,... 25-1. (10) 



One can see that the 25 gapless planes (10^ exist only at 
nonzero J^ , and at J^, = they split off the boundary 
Jr + J^ = to the ferromagnetic phase. 



is for us just a "white spot" which cannot be treated 
within the present approach; to construct an effective 
description for the (S) phase, one has to employ different 
techniques. At finite J^ the gapless line J^ = 2 J^ — Jr 
starts to split off the (BF) boundary J^ = Jr, and the 
(FS) boundary becomes first order, as one sees from the 
behaviour of e-(g) (cf. (||)). The coupling constants 
gA,B,c diverge at the boundaries to the S phase, which 
suggests destruction of any (quasi)-long-range order, and 
thus one may expect that the gapless line terminates at 
the (BS) boundary, though it may in principle continue 
as a first-order transition line. 

It is worthwhile to look at the particular case J^ = Jr. 
According to Ref. |Q, the gapless hue J^ = Jr in this 
case also separates two phases with different string or- 
der, which implies that the lower portion of the B phase 
belongs to the H2 class. It is also known that in this case 
the gapless line continues at larger Jl as the first-order 
line (recall that J^ = J]) = Jr corresponds to the uni- 
form spin chain with next-nearest neighbour interaction) . 
Thus, it becomes clear that additional phase bound- 
aries should exist somewhere inside the spiral "phase," to 
achieve a proper separation of HI- and H2-type phases 
(see Fig. 0(b,c)). This could be an interesting topic for 
the future work. 



IV. SUMMARY 



III. DISCUSSION 



Let us concentrate on the case 5 = ^ as being the 
most important one. For S — -^j a, sketch of the resulting 
phase diagram is presented in Fig. g in a form of two- 
dimensional slices through the phase space at three fixed 
values of J^ {Jr is considered to be positive). 

At J^ = there are no other gapless lines except 
the boundaries of the ferromagnetic phase, and there is 
no spiral phase. The coupling constants gA^ gs diverge 
at the (AB) boundary J^ — 2Jl, which indicates that 
this classical phase boundary gets destroyed by quantum 
fluctuations. On the other hand, all the couphng con- 
stants remain finite at the (BC) and (AC) boundaries, 
but they undergo a jump when crossing the boundaries, 
which suggests a first-order transition. This is in agree- 
ment with the numerical ||ll|,|l^ and bosonization [ [ll|Jl4| 
studies, showing the presence of a first-order transition 
with Jd — Jr/2 being the asymptote for the transition 
line at Ji — !■ oo. There is also a "mirror" transition line 
Jl = Jr/2 due to the Jo ^^ Jl symmetry. According to 
the classification of Ref. [Q , those two first-order tran- 
sition lines separate two topologically different Haldane- 
type phases with Ceven 7^ (phases A,B) and Oodd 7^ 
(phase C); below we refer to those two phases as HI and 
H2 (see Fig. 1(a)). 

At finite J^, the spiral phase (S) appears classically 
in a finite region of the phase diagram. The S region 



We have studied the phase diagram of the generalized 
ladder model with unequal diagonal couplings Jd, J^ 
within the framework of the nonlinear sigma model. We 
show that the phase diagram has a rich structure includ- 
ing several first- and second-order transition boundaries. 
There exist 25 gapless phase boundaries which split off 
the boundary to the ferromagnetic phase at Jd y^ J'^- 
We consider the case 5 = ^ in more detail and show that 
the gapless plane is an extension of one of the transition 
lines discussed in Ref. ||lj] which separate Haldane-type 
phases with different topological order parameter. Still, 
several features of the phase diagram remain unclear and 
require further study. 
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FIG. 1. A schematic representation of the generalized 
ladder model (nl). A, B, C denote different commensurate 
classical ground state configurations. 
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FIG. 2. A sketch of the phase diagram of the model (|l|), 
shown are the slices in the three-dimensional phase space 
at fixed J^ = Jd — J'o'- (a) Jd = 0; (b) Jd ~ Jr\ (c) 
J^ = 3Jr/2. Thick solid and dashed lines denote the sec- 
ond- and first-order transition boundaries, respectively. Thin 
dashed lines indicate crossover between different classical con- 
figurations; the coupling constant diverges at those lines. A, 
B, C are classical configurations shown in Fig. hi F denotes 
the fully saturated ferromagnetic phase, and S stands for the 
spiral "phase" inside which our approach is not valid. HI and 
H2 denote topologically different Haldane-type phases with 
Ceven 7^ and Oodd 7^ accordiug to the classification of Ref. 
0- 



